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Abstract 


We analyze Kl lepton pair decays of Kl —> l + l ~7 and Kl —> l + l~l ,+ l'~ 
( l , V = e, p) within the framework of the light-front QCD approach (LFQA). 
With the Kl —> 7*7* form factors evaluated in a model with the LFQA, 
we calculate the decay branching ratios and find out that our results are all 
consistent with the experimental data. In addition, we study Kl —> l + l~ 
decays. We point out that our prediction on Kl —> e + e“ is about 20% 
smaller than that in the ChPT. We also discuss whether one could extract the 
short-distance physics from Kl —» P + P~ ■ 
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I. INTRODUCTION 


The study of kaon decays has played a pivotal role in formulating the standard model 
of electroweak interactions |]. In particular, the rare decay of K L —> /i + /U was used to 
constrain the flavor changing neutral current [§j as well as the top quark mass 0 . However, 
there are ambiguities in extracting the short-distance contribution since the long-distance 
contribution dominated by the two-photon intermediate state is not well known because its 
dispersive part cannot be calculated in a reliable way | 4 ]-[ 7 j] . To have a better understanding 
of this dispersive part, it is important to study the lepton pair decays of the Kj j meson such 
as Kl —> l + l ~7 and l + l~l + l~ (/ = e,/i) since they can provide us with information on the 
structure of the K L —> 7*7* vertex 0 - 0 . On the other hand, since these lepton pair decays 
are dominated by the long-distance physics, they can also be served as a testing ground 
for theoretical techniques such as chiral Lagrangian or other non-perturbative methods that 
seek to account for the low-energy behavior of QCD. 

Recently, several new measurements of the decay branching ratios of Kl —> P + P~li 
K l —> e + e _ e + e _ , and K L —> e + e~/r + /r~ have been reported [§-[ 11 ]. These decays proceed 
entirely through the K 7*7* vertex and provide the best opportunity for the study of its 
form factor. In Ref. [ 0 ], since the assumption of neglecting the momentum dependence 
for the form factor was adopted, the results for the decays are only valid for those with 
only the electron-positron pair. I11 Ref. 0 ], the decays were studied at the order p 


m 


Chiral Perturbation Theory (ChPT). However, all the results in Ref. pi are smaller than 


the current experimental values. In this work, we consider another non-perturbative method 
in the LFQA to analyze the K 7*7* form factor. As is well known 0 ], the LFQA allows an 
exact separation in momentum space between the center-of-mass motion and intrinsic wave 
functions. A consistent treatment of quark spins and the center-of-mass motion can also be 
carried out. It has been successfully applied to calculate various form factors [[To| [T% . 


The paper is organized as follows. In Sec. II, we derive the theoretical formalism for the 
decay constant and the A"q*7* vertex and use these formalism in the LFQA to extract the 
decay constant and the form factor. I11 Sec. Ill, we fix the parameters appearing in the wave 
functions and calculate the form factors and branching ratios. Finally, conclusions are given 
in Sec. IV. 


II. FRAMEWORK 


We start with the K meson decay constant / A , defined by 


{V\A»\K(P)) = if K P\ 


( 2 . 1 ) 


where 2U = 117^75.8 is the axial vector current. Assuming a constant vertex function A A 
00 which is related to the us bound state of the kaon. Then the quark-meson diagram, 
depicted in Fig. 1 (a), yields 


(0|V|A-(P)} = ~/N c J 


d A Pi A m\ + m s ) i(tfi + rn u ) 

——A A -Tr 

(2 7TJ 4 [ P 2 — m s + * e Pi ~ m u + le 


( 2 . 2 ) 
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where m U:S are the masses of u and s quark, respectively, and N c is the number of colors. We 
consider the poles in denominators in terms of the LF coordinates (p~,p + ,p_ l) and perform 
the integration over the LF “energy” pf in Eq. (| 2 . 6 |). The result is 


( 0 \A»\K{P)) 



[d Pi] __ / j/i I 

Pi P2 P~ - Plan ~ Plan 1 =Pl ° n 


( 2 . 3 ) 


where 


[d 3 Pi] = 


dpi d 2 p 


1 L 


2(2vr) 


Pion 


ml + Pj± 

pT 


Ii = Tr[ 7 5 (^ 2 + m s ) + m u )\. 


( 2 . 4 ) 


For K l —> 7*7*, with the assumption of CP conservation the amplitude is given by 


A(K l -*■ 7*(gi, ei) 7*(?2, £2)) = iF(ql, ql) e pvpa e 2 q[ q%, ( 2 . 5 ) 

where the form factor of F(qf, q f) in Eq. (p. 5 |) is a symmetric function under the interchange 
of ql and q%. In our model, by using the same procedure as above, from the quark-meson 
diagram depicted in Fig. 2 we get 


A(K l - 7* 7*) = - j 
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+ (ei «-»■ e 2 ) 2, (2.6) 


where p 2 = Pi — 9i, P3 = p\ — P, and Cw is the effective contribution to the inclusive s —> d 7* 
decay. After integrating over p±, we obtain 
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( 2 . 7 ) 


where q 2 = P — q\ and 


I 2 = Tr[7 5 (^ 3 + m s ) fa{i> 2 + m s ) /i(j)i + m d )\. 


( 2 . 8 ) 


We note that we do not expect that the absolute decay widths of K L —► Z + Z 7 and —*• 77 

However, we can estimate 


calculated from Eq. (| 2 . 7 |) can fit to the experimental values 
the relative form factors of these leptonic decays versus the two-photon decay, and compare 
the branching ratios with the experimental ones. Recent works on both short-distance (SD) 
and long-distance (LD) contributions to s 
As described in Ref. 


dy* can be found in Ref. [ET 


22 


the vertex function A k l and the denominators in Eq. (| 2 . 7 | ) 
correspond to the Kl meson bound state. In the LFQA, the internal structure of the 
meson bound state ||TT|,|T 8 |,| 23|1 consists of 0 , which describes the momentum distribution 
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s s 

of the constituents in the bound state, and P A (| 2 , which creates a state of definite spin 

(S,S Z ) out of LF hclicity (Ai, A 2 ) eigenstates and is related to the Melosh transformation 

||24|| . A convenient approach relating these two parts is shown in Ref. [^]. The interaction 

Hamiltonian is assumed to be Hi = if PxTysT*!? where T is the quark held and $ is the 

s s 

meson held containing 0 and R x ’ x . When considering the normalization of the meson state 
depicted in Fig. 1 (b) in the LFQA, we obtain 

<M(P',S',S')|P/ Hj\M(P,S,S z )) =2(2n) 3 6 3 (P'-P)S S s'S Sz s' z 

x f [d 3 Pi]0 2 PS 2 Pf;gTr 


7s- 


02 + m 2 01 + mi 

-75' 


P 2 


Pi 


If we normalize the meson state and the momentum distribution function 0 as 117 
<M(P', S', S')|P, H:\M(P, S, S z )) = 2(2vr) 3 P+5 3 (P' - P)S SS 'Ss z S' z , 

and 

f Ppi 1 2 _ 

J 2 ( 2 tt) 3 P +^ 1 

respectively, where pi and p 2 are the on-mass-shell momenta, we have that 


R 


s,s z 

Al,A2 


Pi P 2 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


2VPlon ' P2on + mi?7l 2 

The wave function and the Melosh transformation of the meson are related to the bound 
state vertex function A m by 

A m 


R Xi% § M ■ 


p _ v ~ _ v ~ -”Al,A 2 - (2-13) 

We note that p \, P 2 and p% in the trace of /i >2 must be on the mass shell for self-consistency. 
After taking the “good ” component p = +, we use the definitions of the LF momentum 
variables (x,x',k±,k' ± ) |T 8 | and take a Lorentz frame where P± — P^ — 0 to have q± = 0 
and k' ± = k±. The decay constant Jk and the form factor F{q\,ql) can be extracted by 
comparing these results with Eqs. ( 0 ) and ( j2.5|) , respectively, be., 

dxd 2 k _l 4>K L (x,k±) 


f K = 2V2JN C 


and 


n.UD = j^\cM) j+ 
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where 


a = m Utd x + m s (l - x ), m u = m d , 

r ± = ^r-[ M k + 4 - 4 ± + 9 ? - 44 - 4 Ml L 4 ], 


(2.16) 


and x is the momentum fraction carried by the spectator antiquark in the initial state. 

In principle, the momentum distribution amplitude (j>{x, k±) can be obtained by solving 
the LF QCD bound state equation [^]. However, before such first-principle solutions are 
available, we shall have to use phenomenological amplitudes. One momentum distribution 
function that has often been used in the literature for mesons is the Gaussian-type, 


4>{x, k_ l) g = M\ 


I clk z 
dx 


exp 


k 2 


2 u 2 


(2.17) 


where J\f = 4 ( 7 t/u ; 2 ) 3 / 4 and k z is of the internal momentum k = (k±, fc 2 ), defined through 


1 — x — 


e-i - k z 


e\ + e-2 


x = 


62 + k~ 
ei + e 2 


with d = \/m‘j + k' 2 . We then have 


Mo — e\ + e 2 , k z — 


xM 0 mi + ki 


2 xMn 


and 


dk z 


eie 2 


dx x(l — x)Mq 

which is the Jacobian of the transformation from (x, k±) to k. 


(2.18) 


(2.19) 


( 2 . 20 ) 


III. NUMERICAL RESULTS 


To examine numerically the form factor derived in Eq. 
parameters appearing in (pM{x,k ±). To fit the meson masses, in Ref. |26 


(12.151), we need to specify the 
= 0.22(0.25) 


m,, 


GeV and m s = 0.45(0.48) GeV are obtained with some interaction potentials, while in Ref. 
27| m v = 0.25 GeV and m s = 0.37 GeV in the invariant meson mass scheme. Here we do 


not consider any potential form and scheme and just use the decay constant fx = 159.8 MeV 


0 . charge radius (t 2 )k = 0.34 frn 2 0 . and the quark masses of m Ui d to constrain the s 
quark mass of m s and the scale parameter of uj in Eq. ( eup - By using m u = md = 250 MeV 
0 , we find that m s = 400 MeV and u> = 0.38 GeV. We note that the lower mass of m s 


should not affect the meson masses once we choose a suitable potential or scheme 
Now, we use the momentum distribution functions (p(x, k±)s to calculate the form factors 
F(qf, q 2 ) in time-like region of 0 < q\ and q 2 < M\ ~ 0.25 GeV 2 . In this low energy region, 
we neglect the momentum dependence of the effective vertex Cw(q 2 ) in Eq. (|2.15|), that is, 
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Cw(q 2 ) — Cw{ o). 


( 3 . 1 ) 


We can use Eqs. (|2.15|) and d3d| ) to get the function f(y) = F(qf, 0)/F(0, 0), where y = 
ql/Mx, and the result for \f(y)\ 2 is shown in Fig. 3. From the figure, we see that our result 
with the assumption of Eq. Q3.1|) agrees well with experimental data [3C-32], especially in 
the lower y region. To get a better fit for a larger y, we may use 


C w (q 2 ) ^ 


Cw{ 0 ) 

( 1 - 4 )’ 

V mr. / 


(3.2) 


As seen from Fig. 3, we find that the fit for n (< 8 ) is better than that for n— 1. In particular, 
a larger value of n is preferred if we disregard the data from E845 at BNL |3l| in Fig. 3. 
The experimental result on Kl —■► /r + /i _ 7 from NA48 at CERN, which is currently being 


analyzed |33|], should help to resolve this matter. To illustrate our results on the lepton pair 


decays, we shall take n — 0 and 3, referring as (I) and (II), respectively. 

The function of f(y) is related to the differential decay rate of K® —> l + l ~7 by 


dBm-. _ dT(K L —»■ 7 ) 


dq\ 


/ a 


T(AA —► 77) dqj qj \3n 


\f(y)\ 2 X 3/2 1,774,0 Gi(qf), 


q 2 


Mh 


(3.3) 


where 


A(ci, 6 , c) — q? -j- -I - (? — 2 (ab H- be -I - eti), 


(3.4) 


and 


GM 2 ) = 1 


4 M? \ 


1/2 


1 + 


2 M? 


(3.5) 


Integrating over q\ in Eq. (p.3|) , we get the branching ratios 

T (K° l - e+e- 7 ) 


B e + e 7 — 


B, 




nn - 77 ) 

_ nn - 

Y{Kl - 77 ) 


= 1.64, 1.65 x lO" 2 , 
= 5.50, 6.20 x KT 4 , 


(3.6) 


for (I) and (II), respectively. These values agree well with the experimental data: B g+g-^ = 
(1.69 ± 0.13) x 10 ~ 2 J28[ and B ^+ P Al - 7 = (6T1 ± 0.31) x 10 -4 ||, where we have used [fl4j] 


r exp (j^o 77 ) = [( 5.92 ± 0.15) x 10” 4 ]T exp (A7 -»■ all). 


(3.7) 


On the other hand, our results are larger than B e + e - 7 = 1.59 x 1CT 2 and B I1 + IX - 1 = 4.09 x 10~ 4 , 
respectively, obtained in Ref. JT 2 |, where the momentum dependence of the form factor was 
neglected, i.e., f(y) = 1. This inconsistency is reasonable because the kinematic factor Gi(q 2 ) 
which leads the contribution at q 2 ~ AM 2 is important, and the electron mass is very small 
so that f(y) — 1 is only valid for the decay with an electron-positron pair. For the muonic 
pair case, since the mass of muon is not small, the effect of the deviation of neglecting the 
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momentum dependence is evident. This situation also occurs in the decays with two lepton 
pairs. 

Next, Eq. (|2.15|) can be also used to calculate the differential decay rates of AT —> 
by 


dr(K L ^i + i-r+r-) _ 2 / ck \ 2 
Y(K l 77 ) dql dql q\q\ V 37 tJ 


F(dl4) 


F( 0 , 0 ) 


A3/2 (r Jh ’ Gl(ql)GM) - (3 ' 8) 


After the integrations over q\ and for (I) and (II) we obtain the branching ratios as 
follows: 


Be+e~e+e~ — 
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r(A'» 


e 1 e e 1 e 


B 


T(Kl -> 77 ) 
T(Jt° —> n + e + e~) 

T(K° l -> 


M^M MM 


r(A£ 77 ) 


= 6.61, 6.74 x 10" 5 , 
= 3.87, 4.37 x 10" 6 , 
= 1.50, 1.73 x 10~ 9 . 


(3.9) 


In Table 1, we summary the experimental and theoretical values of the decay branching 
ratios for the K L lepton pair modes. The results of Ref. [jl2] correspond a point-like form 
factor, while those in Ref. |13) are calculated at 0(p d ) in the ChPT. 

Table 1: Summary of the lepton pair decays of AT- 


Br 

PDG j28 

new data 

(I) 

(II) 

Ref. JT3 

Ref. T| 

10 2 x H e + e - 7 

1.69 ±0.09 


1.64 

1.65 

1.59 

1.60 ±0.15 

10 4 x H /t+/J - 7 

5.49 ± 0.49 

6.11 ±0.31 @ 

5.50 

6.20 

4.09 

4.01 ±0.57 

10 5 x B e+e+e - e - 

6.93 ±0.20 

6.28 ±0.65 H 
6.20 ±0.69 10 

6.61 

1 

6.74 

5.89 

6.50 

10 e x B^+^-e+e- 

^• y -4.0 

4.43 ±0.84 11 

3.87 

4.37 

1.42 

2.20 ±0.25 

10 9 x B II+IJ -, J + I1 - 



1.50 

1.73 

0.946 

1.30 ±0.15 


From Table 1, we may also combine the experimental values by assuming that they are 
uncorrelated and we find that 


^, v , = (w±oa)xr‘, 

0"he + ,-e + ,- = (6.83 ± 0.19) xio- 5 , 

= (4.44±°;«) x 1(T«. (3.10) 


It is interesting to see that our results for AT —> l + l ~7 are larger than those in Refs. |]T2|,[l3 
and agree very well with the experimental data. Furthermore, as shown in Eq. ( |3.9| ), those for 
AT — * ► e + e _ e + e _ and AT —*► fi + p~e + e~ also agree with the combined experimental values 
in Eq. Here, we do not consider the interference effect B0] from the identical 

leptons in the final state. The reasons are given in the following. When we use the non- 
point-likc form factor, this effect is about 0.5% in the e + e“e + e~ mode [l3j], which is beyond 
experimental access. For the p + (i~ /i + fi~ mode, the relative size of the interference effect is 
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larger, but it is outside the scope of future experiments because the total branching ratio is 
predicted to be about 8 x 10 -13 . 

We now use the form factor F(qf,g|) to calculate the decays of Kl —■> l + l~. The decay 
branching ratios of the modes can be generally decomposed in the following way 


B l+l - = 


r (k l i+i~) 
T(K l - 77) 


= llm Ai\ 2 + I Re Ai 


(3.11) 


where Irn Ai denotes the absorptive contribution and Re Ai the dispersive one. The former 
can be determined in a model-independent form of 


llm Ai\ = 


a 2 M 2 
2 M* KL fr 


In 


l - A 
l + A 


(3.12) 


where /3f = 1 — 4M 2 . The latter, however, can be rewritten as the sum of SD and LD 

contributions, 


Re Ai = Re A, 


l SD 


Re A, 


l LD- 


(3.13) 


In the standard model, the SD part has been identified as the weak contribution repre¬ 
sented by one-loop IT- box and Z-exchange diagrams |3| j35| , |36[| , while the LD one is related 
to F(q 2 , ql) by 


. l2 2 a 2 M?/3i . . 2 m 2 

I Re x4/ld| = n 2 ]\,f 2 —|Re TZi(M K l )\ , 


(3.14) 


where 37 


Pi(P 2 ) = 


2 i 


d A q 


[P 2 q 2 - (P ■ qf] F(q 2 ,(P-q) 2 ) 


n 2 M\ J g 2 (P - q) 2 [(q - Pl ) 2 - Mf] F(0 , 0) 

In general, an once-subtracted dispersion relation can be written for Re 1Z as 

p2 poo 


(3.15) 


Re 7 Zi(P ) = Re 72./(0) + 


dp ,2 Im Ki(P‘ 


12\ 


7r J o 


f p/2 _ p2^p/2 ’ 


(3.16) 


where Re72-i(0) can be obtained by applying Eq. (|2.15|) in the soft limit of P —> 0. 

For the K L —> e + e _ decay, with n — 0 and 3 of (I) and (II) in Eq. ( |f2l) we find that 


|Re Ae ld| 2 = 5.60, 6.52 x 10” 9 , 

respectively. Since the SD part of ReZUsD can be neglected, we get 

Bi+.- = 1.09 x 10" 8 , 


(3.17) 


B 


e+ e 
II 


= 1.18 x 10' 


(3.18) 


where we have used |Im*4. e | 2 = 5.32 x 10 9 . In terms of the total decay branching ratio 
F> e + e - = T (Kl —> e + e~)/Y{K l —> all), the numbers in Eq. (|3.18|) are about 6.5 and 
7.0 x 10- 12 , respectively. Both results in Eq. (|3.18|) are consistent with the experimental 





























value of £>g+g- = (lAlJ; 9 ) x 10~ 8 measured by E871 at BNL [39|, but they are lower than 
the value of (1.52 ± 0.09) x 10“ 8 [ B e + e - = (9.0 ± 0.5) x 10” 12 ] given by the calculation in 
Ref. |J with the ChPT. It is interesting to note that B e + e - slowly increases as n and reaches 
1.22 x 10~ 8 for n = 10. Clearly, our prediction is about 20% smaller than that in the ChPT 

i- 


For the AT 


/i + /i decay, by subtracting between the value of |Im^4 M | 2 = 1.20 x 10 


-5 


from the experimental data of B e p + p - = (1.21 ± 0.04) x 10 5 |28|.|4(J, we obtain that 

|Re A^ 2 < 7.2 x 10~ 7 (90% C.L.). 

In the standard model, we have that ( 7|,[4T[ 


|Re A fl sD\ 2 Bi< L ^'y'y — 0.9 x 10 9 (1.2 — p) 



3.1 

' ivy' 

170 GeV 


_0.040_ 


(3.19) 


(3.20) 


where p = p(l — A 2 /2). Using the parameters of m t {rrit) = 166 GeV, \V c b\ = 0.041 and 
p ~ 0.224 |J6],(42|, from Eqs. (|T7|) and (|3.20|) we get 


Re A^sd — —1-22 x 10 


-3 


(3.21) 


which is larger than the limit in Eq. (|3.19 ). It is clear that the value of Re*4.^LE> has to be 
either very small for the same sign as Re^gD or the same order but the opposite sign. 

For the case of (I), from Eq. (|3.16|) we find 


Re A 7 LD = -1.11 x 10 


-3 


(3.22) 


which is very close to the SD value in Eq. (|3.21| ) and clearly ruled out if the absolute sign 
in Eqs. (|3.21|) and (|3.22|) are the same. However, if the relative sign is opposite, the limit in 
Eq. (|3.19| ) can be satisfied for certain values of p. From Eqs. ( |3.20| ), ( |3.19| ), and (|3.22| ), by 
taking fh t (rn t ) = 166 GeV and |V%| = 0.041 we extract that 


p > —0.37 or p > —0.38 (90%C.A.). 


(3.23) 


We note that the limit in Eq. ( |3.23|) is close to that in Eq. (41) of Ref. 0. This result is 
not surprising. If we fit F(qf, g|) in Eq. (|2.15| ) with Eq. (14) of Ref. [[7| given by 


f^Ul) = 


F(ql,q 2 2 ) 


— 1 T OL 


q\ 


+ 




ql — m 2 q 2 -rn 2 p/ 


+ (3 


q\q\ 


{ql - m ‘ 2 p ){q 2 - mj) 


F{ 0,0) 

we hnd that a ~ —0.585 and f3 ~ 0.191 and thus 

1 + 2a + /? = 2.16 x 10~ 2 ~ 0 , 

which satisfies the bound of Eq. (35) in Ref. [|7[. Similarly, for (II) we obtain 


Re A “ LD = -1.38 x 10 


-4 


(3.24) 


(3.25) 


(3.26) 


It is very interesting to see that the value in Eq. ( |3.26| ) is much smaller than ReUL^s^ in 
Eq. (|3.21|) , which is exactly the case discussed in Ref. ||. From Eq. ( p.26|) , with the same 
parameters as (I), we hnd that 
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p> 0.63, 0.41 or p > 0.65, 0.42 (90 %C.L.) 


(3.27) 


for the same and opposite signs between ReA^sD and Re A ! /t r Lr) , respectively. We note that 
the limits in Eq. (|3.27|) do not agree with the recent global fitted value of p — 0.224 ± 0.038 
|36| , |42[] , which may not be unexpected since (i) we have not included various possible ranges 
of rh t (mt), | V c fj |, and quark masses in the calculation and (ii) we still need to fix n in Eq. 
( p.2|) and modify the form of CV(g 2 ) |43|]. However, the important message here is that the 
LD dispersive contribution in K Iy —> p, + p~ is calculable in the LFQA. From our preliminary 
results, it seems that R gA^ld is indeed small as anticipated many years ago in Ref. j[§. 
Moreover, our approach here provides another useful tool for the decays beside the ChPT. 


IV. CONCLUSIONS 


In this work, we have studied the Kl lepton pair decays of Kl —> l + l ~7 and Kl —> 
in the light-front QCD framework. In our calculations, we have adopted the 
Gaussian-type wave function and assumed the form of the effective vertex Cw{(f) in Eq. 
( f). 2 |) to account for the momentum dependences in the low energy region. We have calculated 
the relative form factors of the leptonic decays vs. the two-photon decay, and have showed 
that our results on the decay branching ratios of K L —> l + l ~7 and e + e~l + l~ (l = e, p) agree 
well with the experimental data. The remarkable agreements indicate that our form for 
Cw(q 2 ) is quite reasonable, but the number of n still needs to be fixed. Furthermore, all our 
predicted values for these decays are larger than those in the ChPT ||T2| , |T3| 
for the modes of p + p ~7 and p + p~e + e~ for which the 0{jp 6 ) ChPT results in Ref. |13| 
ruled out by the new experimental data |S].[I]J. On the other hand, for Kj J —> e + e _ , we 
have found that B e + e - is between 1.09 and 1.22 x 10 ~ 8 for n = (0,10), which are lower 
than (1.52 ± 0.09) x 10 -8 in the ChPT ||. For Kj j —> p + p~ 1 we have demonstrated that 
the long-distance dispersive contribution is possibly small. However, to get a meaningful 
constraint on the CKM parameters, further theoretical studies [fhj as well as more precise 


in particular 


are 


experimental data such as those from NA48 at CERN B3] on the spectra of the pair decays 
are needed. Finally, we remark that our approach cannot calculate the absolute decay widths 
of Kl —> / + Z - 7 and Kl —■> 77 . 
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FIGURE CAPTIONS 


Fig. 1 Feynman diagrams for the meson (a) decay constant and (b) normalization. 


Fig. 2 Feynman triangle diagrams with (a) and (b) corresponding to the LF valence 
configuration. Empty circles indicate LF wave functions. 

Fig. 3 The y-dependent behavior of \f(y)\ 2 , where the lines from bottom to top correspond¬ 
ing to n = 0,1, • • •, 10 are obtained by this work with fx = 159.8MeV and m s = 400MeV 
and the experimental data are taken from E799 at FNAL | 30|1 , E845 at BNL [31], and NA31 
at CERN |]32|| , respectively. 
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